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Solving for Appendix A - Java Code 

Code for Trajectory Using Law of Sines and Law of Cosines and Pi Space Examples 

Note: This is an old snapshot of the code. The latest code can be viewed and downloaded at https://github.com/pispacesw 

package pispace.core; 

import java.util.ArrayList; 
import java.util.HashMap; 

import java.math.BigDecimal; 
import java.math.Biglnteger; 
import java.math.MathContext; 

import org.apache.log4j.Logger; 

import pispace.core.trajectory.PiSpaceCriteria; 

import pispace.core.trajectory.PiSpaceObjectTrajectorylnfo; 

/* 

* 
* 

* Core class for the Pi-Space Theory. 

* 

* @author Martin Brady 

* 
* 

* This class is a core Pi-Space class with implementations of its 

* important functions. The Pi-Space versions 

* are derived in the Pi-Space Physics Theory in the Advanced 

* Formulas section which contains the complete list. 

* 

* Currently, this class is a work in progress and I will add 

* the functions over time. 

* 
* 

* TERMS AND CONDITIONS: PLEASE READ 

* 

* Although this code may be freely available to use for test applications, 

* this code is NOT free to use in commercial applications and will require 

* payment to the author. The code is based on the Pi-Space Physics Theory 

* which is (Copyright) Martin Brady. 
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public class PiSpaceFormulas { 

static Logger logger = Logger.getLogger(pispace.core.PiSpaceFormulas.class); 

public final double SPEED_OF_LIGHT_KILOMETRES_PER_SECOND = 300000.0 
public final double SPEED_OF„LIGHT_METRES_PER_SECOND = 299792458.0; 
public final double SPEED_OF_LIGHT_MILES_PER_SECOND = 186000.0; 
public final double SPEED_OF_LIGHT_FEET_PER_SECOND = 983571056.0; 

public static final int UNIT_METRES_PER_HOUR = 1; 
public static final int UNIT_MILES_PER„HOUR = 2; 
public static final int UNIT_FEET_PER_HOUR = 3; 
public static final int UNIT_FEET_PER_SECOND = 4; 

public static final double UGCJV1ILES = 3.439E-11; 
public static final double UGCJV1ETRES = 6.67300E-11; 

public int unit; 

/* 

* 

* constructor 

* 

*/ 

public PiSpaceFormulas(int unit) { 
setUnit(unit); 

1 



public static Logger getLogger() { 
return logger; 

} 

public static void setLogger(Logger logger) { 
PiSpaceFormulas.logger = logger; 

} 



/* 

* 

* Metres or Miles? 



Pi-Space Theory, Martin Brady 



Page 3 



*/ 

public void setUnit(int unit) { 
this.unit = unit; 

} 

public int getUnit() { 
return unit; 

} 

/* 
* 

* Hours to seconds 

* 

*/ 

public double getHoursToSeconds(double hours) { 
// return hours; 
return hours / 3600.0; 

} 

/* 

* 

* Seconds to hours 

* 

*/ 

public double getSecondsToHours(double seconds) { 
// return seconds; 
return seconds * 3600.0; 

} 

/* 

* 

* Diameter to energy form 

* 

*/ 

public double getCSquared() { 

if (getUnit() == UNIT_METRES_PER_HOUR) { 

return SPEED_OF_LIGHT_METRESJ J ER_SECOND 

* SPEED_OF_LIGHT_METRES_PER_SECOND 
} else if (getUnitO == UNIT_FEET_PER_SECOND) { 

return SPEED_OF_LIGHT_FEET_PER_SECOND 

* SPEED„OF„LIGHT_FEET_PER_SECOND; 

} else { 

return SPEED_OF_LIGHT„MILES_PER_SECOND 

* SPEED_OF_LIGHT_MILES_PER_SECOND; 

} 
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/* 

* 

* Speed of Light in chosen units 

* 

*/ 

public double getC() { 

if (getUnit() == UNITJViETRESJ'ERJIOUR) { 

return SPEED_OF_LIGHTJVIETRESJ , ER_SECOND; 
} else if (getUnitO == UNIT_FEET_PER_SECOND) { 

return SPEED_OF_LIGHT_FEET_PER_SECOND; 

} else { 

return SPEED_OF_LIGHT_MILES_PER_SECOND; 

} 

} 

/* 

* 

* Map from Pi-Space energy to traditional Newtonian energy 

* 

*/ 

public double getPiSpaceEnergyToNewton(double energy) { 
return energy * getCSquared(); 

} 

/* 

* 

* Convert a Velocity into an energy format 

* 

*/ 

public double getPiSpaceVelocityToEnergy(double velocityPercentLight) { 

return ((velocityPercentLight * velocityPercentLight) / (getCSquared())); 

} 



/* 

* 
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* Pass in a PiSpaceObjectTrajectorylnfo object and this method will calculate the next 

* position by returning a modified copy of the object. 

* 

* position xl, yl cog xl,yl velocity of object kmps angle movement wrt to 

* gravity field acceleration KPMS (e.g. Gravity Field) 

* 

* This algorithm is based on the Orbit document in the Pi-Space Physics 

* Theory 

* 

* Note: The smaller the value of time t, the more accurate the trajectory. 

* 

*/ 

public PiSpaceObjectTrajectorylnfo getNewtonianPiSpaceTrajectoryNextPosition( 
PiSpaceObjectTrajectorylnfo vals 

){ 

getLogger() 

.debugC ITERATION START 

"); 

PiSpaceObjectTrajectorylnfo h = new PiSpaceObjectTrajectoryInfo(); 

double time = vals.getMilliseconds() / 1000.0; 

// Calculate COG distance using COG xl,yl and satellite xl,yl 

double distanceFromCOG = 

getDistanceBetweenTwoPoints( 
vals.getXl(), 
vals.getYl(), 
vals.getXlCOG(), 
vals.getYlCOG()); 

// 

// store values for next iteration 

// 

h.setAngleWRTGravity(vals.getAngleWRTGravity()); 

h.setMilliseconds(vals.getMilliseconds()); 

h.setXlCOG(vals.getXlCOG()); 

h.setYlCOG(vals.getYlCOG()); 

h.setRadius(vals.getRadius()); 
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h.setMassl(vals.getMassl()); 

h.setMass2(vals.getMass2()); 

h.setCriterias(vals.getCriterias()); 

h.setCurrentVelocity(vals.getCurrentVelocityO); 

h.setlnitialAngleWRTGravity(vals.getlnitialAngleWRTGravityO); 

getLogger().debug("Object initial velocity is " + vals.getlnitialObjectVelocityO); 
getLogger().debug("Milliseconds is " + vals.getMillisecondsO); 
getLogger().debug("Time milliseconds/1000 is " + time); 

getLogger().debug("AngleAxesOffset from field is " + vals.getAngleAxesOffset()); 
getLogger().debug("DistanceFromCOG is " + distanceFromCOG); 
getLogger().debug("Radius of planet is " + vals.getRadius()); 
getLogger().debug("AngleWRTGravityField is " + vals.getAngleWRTGravityO); 
getLogger().debug("Initial AngleWRTGravityField is " + vals.getlnitialAngleWRTGravityO); 
getLogger().debug("Mass of planet is " + vals.getMass2()); 
getLogger().debug("Criterias are " + vals.getCriterias().size()); 
getLogger().debug("Current Field Velocity is " + vals.getCurrentVelocityO); 
getLogger().debug(""); 



// 

// velocity based distance for field movement 

// 

double vt = 0; 

double fieldVelocity = 

(vals.getlnitialObjectVelocityO 

* 

Math.cos(convertDegreesToRadians(vals.getInitialAngleWRTGravity()))*(-1.0)); 
h.setlnitialConstantFieldVelocity(fieldVelocity); 

getLogger().debug("Initial constant field velocity is " + h.getlnitialConstantFieldVelocityO); 

if (vals.isUsePiSpaceFormulasO) { 

vt = this.getPiSpaceDistanceUseNewtonianUnits(fieldVelocity, 0, time, true); 

} else { 

vt = this.getNewtonianDistance(fieldVelocity, 0, time); 

} 

getLogger().debug("V*T field distance is " + vt); 
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// velocity based distance for particle movement, does not change 

double particlevt = 0; 

double particle Velocity = 

(vals.getlnitialObjectVelocityO 

* 

Math.sin(convertDegreesToRadians(vals.getInitial Angle WRTGravityO))); 
h.setInitialConstantParticleVelocity(particle Velocity); 

getLogger().debug("Initial constant particle velocity is " + h.getlnitialConstantParticleVelocityO); 

if (vals.isUsePiSpaceFormulasO) { 

particlevt = this.getPiSpaceDistanceUseNewtonianUnits(particle Velocity, 0, time, true); 

} else { 

particlevt = this.getNewtonianDistance(particle Velocity, 0, time); 

} 

getLogger().debug("V*T particle distance is " + particlevt); 
getLogger().debug(" "); 

// determine gravity at this point 

double grav = getNewtonianGravity(vals.getMassl(), vals.getMass2(), distanceFromCOG); 
getLogger().debug("Derived gravity for this position is "+grav); 

//Y velocity relative to gravity field e.g. straight up 

//We need to keep track of this from the previous moment 

double velocitylnYDirection = this.getNewtonianFinalVelocity( 

vals.getFieldVelocityO, grav, vals.getFieldDistance()); 

getLogger().debug( 

"Field velocity for next step is "+ 
velocity In YDirection+ 

" was "+ 

vals.getFieldVelocity()+ 

" for distance "+vals.getFieldDistance()); 



Pi-Space Theory, Martin Brady 



Page 8 



10/26/2013 



h.setFieldVelocity(velocitylnYDirection); 

// just focus on distance due to field acceleration and time 

double accDist = 0; 

if (vals.isUsePiSpaceFormulasO) { 

accDist = this.getPiSpaceDistanceUseNewtonianUnits(velocityInYDirection, grav, time, true); 

} else { 

accDist = this.getNewtonianDistance(velocityInYDirection, grav, time); 

} 

h.setFieldDistance(accDist); 

getLogger().debug("Acceleration field distance is " + accDist); 

// 

// path of least time... 

// do we continue to move in the direction we are travelling in? or does gravity take over? 
// 

getLogger().debug(" "); 

double pathOfLeastTimeAngle = this.getPiSpacePathOfLeastTime( 
vt, accDist, vals.getAngleWRTGravityO); 

vals.setAngleWRTGravity(pathOfLeastTimeAngle); 

// **** 

// **** st e p x : Calculate distance u and new angle WRT gravity Alpha 

// **** 

getLogger().debug(""); 
getLogger() 

.debug("Step 1: Calculate distance u and new angle WRT gravity Alpha"); 
// First calculate u = new velocity per second 

// Rule of thumb: Extract the field first, then apply the angle WRT to 

// gravity, produces 180-angleWRTGravity 

// Reason to do it this way, we get acceleration in the direction of 

// movement towards COG 

// while using the Law of the Cosines 
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double interiorAngle = 180 - vals.getAngleWRTGravityO; 

getLogger().debug("Interior angle is "+interiorAngle); 
getLogger().debug(" "); 

double fieldDistance = vt - accDist; 

double u = getLawOfCosinesDistance(particlevt, fieldDistance, interiorAngle); 
getLogger().debug( 

"Distance u which is distance between particle lengths " + particlevt 

+ " and field length " + (fieldDistance) + " having interior angle " + interiorAngle 

+ " is " + u); 

getLogger().debug("u is "+u); 

getLogger().debug(""); 

//for now only consider y velocity 

double thevelocitylnYDirection = this.getNewtonianFinalVelocity( 

h.getFieldVelocityO, grav, h.getFieldDistance()); 

getLogger().debug("The current opposing field velocity is "+thevelocityInYDirection); 

h.setCurrentVelocity(h.getlnitialConstantFieldVelocityO-thevelocitylnYDirection); 
getLogger().debug("Current field velocity is (part which changes) "+h.getCurrentVelocity()); 
getLogger().debug(" "); 

// 

//Calculate the Alpha angle 

// 

double pureAlphaAngle = getLawOfSinesAngle(u, interiorAngle, particlevt); 
getLogger().debug("Alpha angle is "+pureAlphaAngle); 
getLogger().debug(" "); 

//initial velocity stays constant 

h.setlnitialObjectVelocity(vals.getlnitialObjectVelocityO); 

// Next calculate the new angleWRTGravity called alpha 

double betaAngle = getLawOfSinesAngle(u, interiorAngle, fieldDistance); 

getLogger().debug("Beta angle is "+betaAngle); 

getLogger().debug(" "); 

// Use accDist, it's better to use than v*t... 

double alphaFromBeta = 180 - betaAngle - interiorAngle; 
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getLogger().debug("Setting new angleWRTGravity (oldAnglelnewAngle) " 
+ vals.getAngleWRTGravityO + "/" + alphaFromBeta); 

// new angleWRTGravity 
h.setAngleWRTGravity(alphaFromBeta); 

// **** 

// **** Step 2: Calculate length r which is the distance to the COG, 

// from this calculate the new length s which is the new distance to COG 

// **** 

getLogger().debug(""); 
getLogger() 

.debugC'Step 2: Calculate length r which is the distance to the COG, from this 
calculate the new length s which is the new distance to COG"); 

double r = distanceFromCOG; 

// 

// Simple collision detection 

// 

if (accDist > r) { 

getLogger() 

.debug("!!!!Gravity pulled object through the COG origin!!!! Collision 

detection."); 

// Right now this is how we can end 
h.setCollisionDetection(true); 

} else { 

h.setCollisionDetection(false); 

} 

// 

// new distances 

// 

double s = getLawOfCosinesDistance(u, r, alphaFromBeta); 
h.setDistanceFromCOG(s); 

getLogger().debug("(r(old distanceFromCOG)ls(new distanceFromCOG)) " 
+ r + "/" + s); 

// **** 

// **** st e p 3. Calculate new angle WRT to axes, new theta 

// **** 
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getLogger().debug(""); 

getLogger().debug("Step 3: Calculate new angle WRT to (field) axes, new theta"); 
double newAngleWRTAxes = getLawOfSinesAngle(s, alphaFromBeta, u); 

getLogger().debug("(angleWRTFieldAxeslnewAngleWRTFieldAxes) " + vals.getAngleAxesOffset() 
+ "/" + newAngleWRTAxes); 

double totalOffset = (vals.getAngleAxesOffset() + newAngleWRTAxes) % 360; 
getLogger().debug("Total axis offset due to particle movement " + totalOffset); 
getLogger().debug(" "); 

h.setAngleAxesOffset(totalOffset); 

// final step 4 

// calculate the new x,y position of the object moving the 
// gravity field 

// we have an angle and a length (polar co-ordindates) so we can covert 
// into 

// a new x,y position. 

double r Angle = convertDegreesToRadians(90 - totalOffset); 

getLogger().debug("X (sl90-totalOffsetlMath.cos) " + s + "/" 

+ (90 - totalOffset) + "/" + Math.cos(r Angle)); 

getLogger().debug("Y (sl90-totalOffsetlMath.sin) " + s + "/" 

+ (90 - totalOffset) + "/" + Math.sin(r Angle)); 
getLogger().debug(" "); 

double newx = s * Math.cos(rAngle) + vals.getXICOGO; 
double newy = s * Math.sin(r Angle) + vals.getYlCOGO; 

getLogger().debug("newx relative to fixed axis " + newx); 

getLogger().debug("newy relative to fixed axis " + newy); 

getLogger().debug(" "); 

h.setXl(newx); 

h.setYl(newy); 

double deltax = newx - vals.getXl(); 
double deltay = newy - vals.getYl(); 

getLogger().debug("deltax " + deltax); 
getLogger().debug("deltay " + deltay); 
getLogger().debug(""); 
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h.setDeltaY(deltax); 
h.setDeltay (deltay ) ; 

h.setTotalDeltaX(h.getDeltaX()+vals.getTotalDeltaX()); 
getLogger().debug("Total Delta X is " + h.getTotalDeltaXO); 
h.setTotalDeltaY(h.getDeltay()+vals.getTotalDeltaY()); 
getLogger().debug("Total Delta Y is " + h.getTotalDeltaY()); 
getLogger().debug(" "); 

double rememberTime = vals.getTotalTime(); 
h.setTotalTime(vals.getTotalTime()+time); 

getLogger().debug("Total time is "+h.getTotalTime()+" ("+rememberTime+"+"+time+")"); 
getLogger().debug(" "); 

h.setDeltaDistance(u) ; 

getLogger().debug("Delta distance from COG is "+h.getDeltaDistance()); 

getLogger().debug("Previous distance from COG "+vals.getTotalDistance()); 

h.setTotalDistance(vals.getTotalDistance()+h.getDeltaDistance()); 
getLogger().debug("Total distance is "+h.getTotalDistance()+" adding delta distance 
"+h.getDeltaDistance()); 

getLogger().debug(""); 

h.setDeltavelocity(h.getCurrentVelocity() - vals.getCurrentVelocityO); 
String status=""; 

if (h.getCurrentVelocityO > 0 && h.getDeltavelocityO > 0.0) { //positive increase 

status="SPEEDS.UP POSITIVE VELOCITY"; 
} else if (h.getCurrentVelocity() > 0 && h.getDeltavelocity() < 0.0) { 

status="SLOWS.DOWN POSITIVE VELOCITY"; 

} else if (h.getCurrentVelocityO < 0 && h.getDeltavelocity() < 0.0) { //negative increase 

status="SPEEDS.UP NEGATIVE VELOCITY"; 
} else if (h.getCurrentVelocity() > 0 && h.getDeltavelocity() > 0.0) { 

status="SLOWS.DOWN NEGATIVE VELOCITY"; 

} else { 

status="STOPPED!"; 

} 

getLogger().debug(" "); 



if (vals.getTotalDistance() != 0) { 
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getLogger().debug("Delta velocity is "+status+" "+h.getDeltavelocity()+"(+speeds.up - 

slows.down) (new:" 

+h.getCurrentVelocity()+"-old:"+vals.getCurrentVelocity ()+")"); 

} 



getLogger().debug(" ITERATION END, CHECK 

CRITERIAS "); 

checkTrajectoryCriterias(vals.getCriterias(), h); 
return h; 

} 



/* 
* 

* check the criterias 

* 

*/ 

public boolean checkTrajectoryCriterias(ArrayList<PiSpaceCriteria> criterias, PiSpaceObjectTrajectorylnfo t) { 

boolean met = false; 

for (PiSpaceCriteria p : criterias) { 
p.perform(t); 
if (p.isCriteriaMet()) { 

met = true; 

p.report(t); 

} 

} 

return met; 

} 



/** 

* @param args 

* 

*/ 

public static void main(String[] args) { 
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@Test 

public void testDistanceCalculationJustGoingUpO { 

//This does not test going up and then coming back down, just going up 

PiSpaceFormulas pse = new PiSpaceFormulas( 

PiSpaceFormulas.UNiT_METRES_PER_HOUR); //I will fix up units later 

getLogger().debug("180 degrees Straight up, use Pi-Space Formulas"); 

double massl = 100; 

double timems = 1000; 

//Planet p = new Planet(5.98E24,6.378E6,"EARTH"); 
//planets.put("EARTH", p); 

double massOfEarth = 5.98E24; 

double radiusOfEarth = 6.378E7; 

double distanceFromCOG = radiusOfEarth; 

double steps = 5; 

double velocity = 100.0; 

double earthGravity = 9.809637070728721; 

//Use the trajectory 

PiSpaceObjectTrajectorylnfo t = new PiSpaceObjectTrajectoryInfo( 
0, distanceFromCOG, //xl,yl position 
0, 0, //x2,y2 cog 
velocity, //velocity 
180, //direction of movement 
timems, //milliseconds 
radiusOfEarth, //radius of planet 
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massl, //mass of object 
massOfEarth, //mass of earth 
true //use Pi-Space formulas 

); 

for (int i = 0; i < steps; i++) { 

getLogger().debug("STEP"+(i+l)); 

t = pse.getNewtonianPiSpaceTrajectoryNextPosition(t); 

} 

getLogger().debug("Total distance summing trajectory is "+t.getTotalDistance()); 

//Use the Newtonian approach, summing it all up in one function 

double vt = pse.getPiSpaceDistanceUseNewtonianUnits( 
velocity, -earthGravity, steps, true); 

getLogger().debug("Total Pi-Space function distance is "+vt); 

vt = pse.getNewtonianDistance( 

velocity, -earthGravity, steps); 

getLogger().debug("Total Newtonian function distance is "+vt); 

double accuracy = 0.1; 

assertEquals( "Passed " , 

t.getTotalDistance(), 
vt, 

accuracy); 



) 



